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ABSTRACT 

Recent advances in observational cosmology are changing the way we view the nature 
of time. In general relativity, the freedom in choosing a time hypersurface has ham- 
pered the implementation of the theory Fortunately, Hamilton- Jacobi theory enables 
one to describe all time hypersurfaces on an equal footing. Using an expansion in 
powers of the spatial curvature, one may solve for the wavefunctional in a semiclas- 
sical approximation. In this way, one may readily compare predictions of various 
inflation models with observations of microwave background anisotropics and galaxy 
clustering. 
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1 Introduction 



In the standard cosmological formulation of Einstein's theory of gravity, one ordi- 
narily solves the field equations for the 4-metric g^. In a more elegant although equiv- 
alent approach, it has been recommended [l]-[6] that one solve the Hamilton- Jacobi 
(HJ) equation for general relativity which governs the evolution of the generating 
functional S. By adopting this method, many advantages have be gained: 

(1) Avoiding Gauge Problems. Neither the lapse N nor the shift Ni appear in 
the HJ equation. The generating functional depends only on the 3-metric 7^ (as well 
any matter fields that may be present). Hence, the structure of the HJ equation is 
conceptually simpler than that of the Einstein field equations. For instance, one is 
able to avoid (temporarily) the embarrassing problem of picking a gauge in general 
relativity. As a corollary, one obtains a deep appreciation for the Nature of Cosmic 
Time [7]. 

(2) Solution of Constraint Equations. One may solve the constraint equations 
of general relativity in a systematic manner, even in a nonlinear setting [4], [8]. The 
momentum constraint is easy to solve using HJ theory: one simply constructs S using 
integrals of the 3-curvature over the entire 3-geometry. The energy constraint may 
be solved by expanding in powers of the 3-curvature ('spatial gradient expansion'). 

(3) Primitive Quantum Theory of Gravity. Solutions of the HJ equation may be 
interpreted as the lowest order contribution to the wavefunctional for an expansion in 
powers of ft, (semiclassical approximation). One may describe some quantum processes 
such as the initial 'ground state' of the Universe [9]-[10] or tunnelling through a 
potential barrier [11]. 

If one accepts that fluctuations for galaxy formation as well as microwave back- 
ground anisotropics were generated during an inflationary epoch, then it is imperative 
that one quantize the gravitational field [9]- [12]. However, quantization of the full 
gravitational field has proven to be elusive. Several possible forms for a viable theory 
have been advanced, although there is no general consensus. String theory is the 
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most popular candidate, and its status has been reviewed by many of the speakers 
in this conference [13]- [14]. The goal here will be more modest. Beginning with the 
HJ equation, I will be content to consider the semiclassical theory of Einstein gravity. 
In this way, one follows in spirit the historical development of the theory of atomic 
spectra. Before the development of the quantum theory in 1926, the semiclassical 
theory of Bohr and Sommerfeld provided a useful although imperfect description of 
various atoms. 

HJ theory has proven to be a particularly powerful tool for the cosmologist. In 
this article, I will review our current understanding of the theory. In addition, I will 
discuss the observational status of three models of cosmological inflation: 

(1) inflation with an exponential potential ('power-law inflation') [15], [16] which 
arises naturally from Induced Gravity [17], [18] or Extended Inflation [19]; 

(2) inflation with a cosine potential ('natural inflation') [20] where the inflaton is a 
pseudo-Goldstone boson; 

(3) inflation with two scalar fields ('double inflation') [21], [17] where there are two 
periods of inflation. 



2 Hamilton-Jacobi Theory for General Relativity 



The Hamilton-Jacobi equation for general relativity is derived using a Hamiltonian 
formulation of gravity. One first writes the line element using the ADM 3+1 split, 



where N and TV; are the lapse and shift functions, respectively, and 7^ is the 3-metric. 
Hilbert's action for gravity interacting with a scalar field becomes 
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The lapse and shift functions are Lagrange multipliers that imply the energy con- 
straint 7i(x) = and the momentum constraint 7ii(x) = 0. 



The object of chief importance is the generating functional 

S = S[)ij{x),4>{x)]. (3) 

For each scalar field configuration <f>(x) on a space-like hypersurface with 3-geometry 
described by the 3-metric 7^(2), the generating functional associates a complex num- 
ber. The generating functional is the 'phase' of the wavefunctional in the semiclassical 
approximation: 

- e lS . (4) 

(The prefactor is neglected here although it has important implications for quantum 
cosmology [24].) The probability functional, V = \^\ 2 , is given by the square of the 
wavefunctional. Replacing the conjugate momenta by functional derivatives of S with 
respect to the fields, 



6'j ij (x) ' S<j){x) 

and substituting into the energy constraint, one obtains the Hamilton- Jacobi equation 
[22], [23], 

+ k 1/2 (^) 2 + ^ 2 ™ 

-\i ll2 R + \i ll2 iHA, = ^ (6) 

which describes how S evolves in superspace. R is the Ricci scalar associated with 
the 3-metric, and V(<f>) is the scalar field potential. In addition, one must also satisfy 
the momentum constraint 

Hi{x) = -2 7 iJfc - — — + - — --7ifc,i + j~r^4>,i = , 7 

which legislates spatial gauge invariance: S is invariant under reparametrizations of 
the spatial coordinates. Since neither the lapse function nor the shift function appears 
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in eq.(6) and eq.(7), the temporal and spatial coordinates are arbitrary: HJ theory 
is covariant. (Units are chosen so that c = 8kG = Ti = 1. The HJ equation for 
Brans-Dicke gravity has been studied in ref.[25].) 

It is quite important that the 3-metric "fij(x) be positive definite otherwise two 
fatal disasters would occur: (1) the principle of microscopic causality would be vio- 
lated; (2) points x and y appearing in the generating functional would no longer be 
space-like, and hence the fields amplitudes <f>(x) and <f>(y) may not be independent 
— - this would mess up the computation of functional derivatives appearing in eq.(6) 
and eq.(7). 



2.1 Partial Reduction of Einstein's Equations 



From a real solution S of the HJ equation and the momentum constraint, one may 
construct solutions to Einstein's equations. Provided that the fields evolve according 
to the definitions of the momenta, 

(<£ - JVV,i) IN = 7" 1 / V, (8) 

(70- " N HJ ~ N j\i) l N = 27" V V (2 7 ^7*/ " HHm) , (9) 

one may verify that Einstein's equation are indeed satisfied. (Note that | denotes 
covariant differentiation with respect to the 3-metric.) As a result, HJ theory allows 
for a partial reduction of Einstein's equations. The lapse and shift appear only in the 
definitions of the momenta for gravity and matter, eq.(8) and eq.(9). They are arbi- 
trary and they reflect the observer's freedom in choosing his space-time coordinates. 
For example, if one chooses time hypersurfaces such that the scalar field is uniform, 
<j) = t, then the lapse is given automatically by eq.(8), 

= (10) 

Typically one sets the shift to zero, in which case the spatial coordinates of one time 
slice are projected orthogonally onto the others. In order to complete the determina- 
tion of the 4-metric as a function of time <f>, one need only integrate eq.(9). 
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2.2 Evolution and Observation of a Gravitational System 



The above discussion demonstrates how HJ theory enables one to split the analysis 
of any gravitational system into two stages: 

(1) Gauge-independent evolution of the system which is governed by the HJ equation 
as well as the momentum constraint equation and 

(2) Observation of the system which is gauge-dependent (i.e., it depends on arbitrary 
choices for N,Ni). 

If the generating functional S is complex, then one is describing intrinsically quan- 
tum processes. The freedom in choosing the lapse and shift then reflect the necessity 
for gauge-fixing. Otherwise, expectation values for physical observables such as the 
two-point correlation function [1], 

< <Hx)<t>(y) >= J [d 7 ][#] <K*)<Ky) l*l 2 , (ii) 

would be infinite. In the spirit of Dirac quantization [23], gauge-fixing is performed 
only after one has determined a solution of the HJ equation (or its quantum analog, 
the Wheeler-DeWitt equation [26]). 



3 Prototype Hamilton-Jacobi Solution 



Until recently, it had been thought that the Hamilton-Jacobi equation for general 
relativity was intractable. As a result, since the pioneering work of DeWitt [26] and 
Misner [27] in the 1960's and 1970's, much effort was spent studying minisuperspace 
models where one effectively truncates superspace: one examines only a finite num- 
ber of variables describing typically a homogeneous universe. For example, in the 
context of quantum cosmology, these models were studied by Hartle, Hawking, Page, 
Vilenkin and others in the 1980's. One hoped to understand some qualitative features 
of semiclassical gravity from these 'toy models'. However, it is generally agreed that 
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such models represent only the initial stages of a much larger program. For example, 
in order to appreciate more fully the nature of time, one must incorporate the role 
of inhomogeneities. After all, a time hypersurface represents the arbitrary manner 
in which one slices a 4-geometry [2]. Many properties of full superspace can be un- 
derstood by employing a series solution method. The prototype solution, the spatial 
gradient expansion, was suggested by John Stewart and myself [5]. Effectively, one is 
able to decompose superspace into an infinite but discrete number of minisuperspaces 
which are tractable. 



3.1 Spatial Gradient Expansion 



As a first step in solving eq.(6) and eq.(7), one expands the generating functional 

S = SW+S&+SW + ... , (12) 

in a series of terms according to the number of spatial gradients that they contain. 
The invariance of the generating functional under spatial coordinate transformations 
suggests a solution of the form, 

<S (0) [7^),<^)] = -2 / ^x^ 2 H[<f>(x)] , (13) 



for the zeroth order term S^°\ The function H = H{<f)) satisfies the separated HJ 
equation (SHJE) of order zero [6], 




^=okr +o^) , (14) 



which is an ordinary differential equation. (The term 'separated' describes the fact 
that the metric variables are absent in eq.(14) — they have been separated out.) 
Note that S (0) cont ains no spatial gradients. The zeroth order term is an excellent 
approximation for universes where the wavelength of any inhomogeneities is much 
larger than the Hubble radius In fact, the notion of long-wavelength field is an 

essential ingredient of any model of cosmic structure formation. 
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3.2 Higher Order Terms 



An important simplification occurs for higher order terms, 

5 (2n) ) f or „ > I- they are 
governed by linear partial differential equations of the inhomogeneous type: 

_ 2 ^ + M «^ + R .* )=0 . (15) 
dcf) 6(f) Sjij 

The remainder term TZ^ 2 ^ depends on some quadratic combination of the previous 
order terms [3]. For example, for n = 1, it is 

= i 7 1/2 7%^- - \i ll2 R , (16) 

whereas for n > 2, the remainder, 7?S 2n \x), is given by 

n-l cq(2p) c c(2n-2p) 

, ™ _1 1 /)5( 2p ) /)5( 2 ™- 2 p) 
+ -i/2 6* (17) 



In order to compute the higher order terms from eq.(15), one introduces a change 
of variables, (7^, 0) — ► ( fij,u): 



d( f> f _ Q -2( N 

where the conformal factor Q = is defined through 

du d(j) d(j) 

The equation for 5^ 2 ™) simplifies considerably: 

65(2") 



8u(x] 



+ 1Z^[f ij (x),u(x)}=() 



(18) 



(19) 



(20) 



Eq.(20) has the form of an infinite dimensional gradient. Before integrating it, I will 
review some elementary results from potential theory 
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3.3 Potential Theory 



The fundamental problem in potential theory is: given a force field g l (uf,) which is 
a function of n variables uj., what is the potential $ = (if it exists) whose 

gradient returns the force field, 



8U: 



(21) 



Not all force fields are derivable from a potential. Provided that the force field satisfies 
the integrability relation, 



= 



d_l_ 



d d 



8ua ' du,; 



(22) 



(i.e., it is curl-free), one may find a solution which is conveniently expressed using a 
line-integral 



$K) = J c J2 dv j ( A v i) 



(23) 



If the two endpoints are fixed, all contours return the same answer. In practice, one 
employs the simplest contour that one can imagine: a line connecting the origin to 
the observation point u^. Using s, < s < 1, to parameterize the contour, the 
line-integral may be rewritten as 



n „i 
j = l ° 



(24) 



3.4 The Nature of Cosmic Time 



In solving eq.(20) for the generating functional S^ 2 ^ of order 2n, one utilizes a line- 
integral in superspace: 

S (2n) = - J d 3 x £ ds u(x) n {2n) [f il (x),su(x)} . (25) 

For simplicity, the contour of integration was chosen to be a straight line in superspace. 
As long as the end points are fixed, the line integral is independent of the contour 
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choice which corresponds to a specific time foliation. This property goes a long way 
in illuminating the nature of time for semiclassical relativity. (Many questions remain 
concerning the role of time in a quantum setting [28]; because the above arguments 
have been quite general, I conjecture that line-integrals in superspace will useful for 
a full quantum description.) 



[2], 



One may verify the integrability of eq.(20) by explicitly computing the commutator 
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8u(x) ' 8u(y) 



5 (2") = 



)u(y) 



8u{x) 



d 



= hHx) Uf n - 2 \x)+^{y) Hf n - 2 \y)\ ^-6 3 (x - y) 



(26) 



which assumes by induction that S^ 2 \ S^, . . . , S^ 2n 2 ^ satisfy eq.(20). The 'inte- 
grability condition' of potential theory, eq.(22), demands that the commutator (26) 
vanish. In the above expression, Hf n ~ 2) is the moment um constraint evaluated using 
the generating functional of order (2n — 2): 



H 



(2*1-2), 



/ ^(2«-2)\ sS (2n - 2) SS {2n ~ 2) 
-2 ljk—c — ]—r + ~ 7 — r-rlki,j + 777 \ <t>,; 



(27) 



We conclude that <S( 2 ™) is indeed integrable provided the term of previous order, 
is invariant under reparametrizations of the spatial coordinates: 7{ 



(2n-2) _ 



0. 



In general, the integrability condition for the Hamilton- Jacobi equation follows from 
the Poisson brackets [29] between the energy densities evaluated at the two spatial 
points x and y: 



{H(x),H(y)} = [^(x)n j (x) + 1 i ^y)H 1 (y)} ^S 3 (x - y) 



(28) 



Typically, <S( 2 ™) is an integral of terms which contain the Ricci tensor and spatial 
derivatives of the scalar field [3]. For n = 1, one determines that 

S^[f tJ (x),u(x)} = J d\f 1 ' 2 [j(u)R + k(u)f ij u iUd ] . (29) 
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R is the Ricci scalar of the conformal 3-metric appearing in eq.(18). The u- 
dependent coefficients j and k are, 

j(u) = j du' + F, k(u) = H{u) n(u) , (30) 

where F is an arbitrary constant. 



3.5 Characteristics of Cosmic Time 



The generalization of the spatial gradient expansion to multiple scalar fields is non- 
trivial [2]. In this case, one employs the method of characteristics for solving the 
linear partial differential equation that appears. For a single scalar field <f> in a HJ 
description, it is more or less obvious to use some function of <f> as the integration pa- 
rameter. In order to facilitate the integration of for multiple fields, I recommend 
using the scale factor, ft = Q((f) a ), which is a specific function of the scalar fields. A 
brief summary follows. 



One first solves the SHJE of order zero, eq.(14), describing two scalar fields to 
find a Hubble function 

H = H^fafcM (31) 

which depends on two homogeneous parameters 4>\ and 4>2- One then makes a change 
of variables 

[4>i{x),(f)2{x),')ij{x)} -> [Q.(x),e(x)Jij(x)\ (32) 

where the new variables are found by computing partial derivatives of H with respect 
to the parameters: 

f f)H\ ~ 1//3 

n«.)= y , (33) 

em= f/|f, (34) 

fij = ^~ 2 (4>a) lij ■ (35) 
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In terms of the new fields, the generating functional of order two can be computed 
explicitly: 



5 (2) = J d 3 x / 1/2 e)R + ku(Sl, e)a,fF + 2k 12 (£l, e)Q. 4 e ;i + k 22 {VL, 
j, kn, k\2 and k 22 , are functions of (fi, e) which are found to be: 

>< n - e >= f^2if(k) +J » (e) ' 



e e,, : e' 



A; 22 (^, e) 
A;i 2 (fi, e) 



/o 

0; 



rffi' n(fi',e) + k (e) , 



(36) 

(37) 

(38) 

(39) 
(40) 



jo = jo(e) and k$ = ka(e) are arbitrary functions of e; n(Ct,e) is defined according to 



n(Q, e) = 



+ 



de Jn \ de Jn 



(41) 



Details are described in ref. [2]. 



4 Quadratic Curvature Approximation 



In order to describe the fluctuations arising during the inflationary epoch, it is nec- 
essary to sum an infinite subset [1], [2], of the terms S( 2n \ In this case, one makes 
an Ansatz which includes all terms which are quadratic in the Ricci tensor R^ of the 
conformal 3-metric fij(x): 

S = + + Q ; (42) 
here the quadratic functional Q is 



Q 



d^xf 1 ' 2 



R S(u, D 2 ) R + R ij f(u, D 2 ) R l} - -R f(u, D 2 ) R 



(43) 



where S(u, D 2 ) and T(u, D 2 ) are differential operators which are also functions of u. 
S and T describe scalar and tensor fluctuations, respectively. D 2 is the Laplacian 
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operator with respect to the conformal 3-metric. Terms which are cubic and higher 
are neglected. 



4.1 Multiple Fields: Quadratic Constant Approximation 

Once again, the case for two fields [2] is more complicated (after which the extension 
to any additional fields is straightforward). One replaces the scalar operator S by a 
matrix operator S a b, a,b = 1,2, which is a function of Q(x) and e(x), eq.(33) and 
eq.(34). The scalar operator S a b is then sandwiched between the vector [R 1 D 2 e] and 
its transpose in the generalization of eq.(43) to the quadratic constant approximation. 



5 Comparison with Large-Angle Microwave Back- 
ground Fluctuations and Galaxy Correlations 



Using H J theory, I will compare the cosmological implications of three inflationary 
models: Model 1 — power-law inflation; Model 2 — natural inflation, and Model 
3 — double inflation. All models will be normalized using large angle microwave 
background anisotropics determined by COBE [30]: 0^(10°) = 30.5 ± 2.7 \xK (68% 
confidence level). 

It is conventional to parametrize the primordial scalar fluctuations arising from 
inflation by ( which is proportional to the metric perturbation on a comoving time 
hypersurface [6]. The power spectra for ( are shown in Figs. (la), (2a), (3a), for 
Models 1, 2, 3, respectively. Both power-law inflation and natural inflation yield 
power spectra V((k) for ( that are pure power-laws: 

/ , \ n-l 

V c (k) = V c (k )[-} . (44) 
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The spectral index for scalar perturbations is denoted by n, and n = 1 describes the 
flat Zel'dovich spectrum. The simplest models arising from inflation are characterized 
by n < 1. The normalization of the spectra differs for the first two models since 
gravitational waves may contribute significantly to COBE's signal for the power-law 
inflation model [16]. The primordial power spectrum for double inflation is not a pure 
power-law. (In Fig. (3a), the primordial fluctuations for inflation with a single field 
having a quadratic potential is also shown.) 

The power spectra Vs(k) for the linear density perturbation 6 = Sp/p at the 
present epoch are shown in Figs. (lb), (2b), (3b). The data points with error bars are 
determined from the clustering of galaxies [31]. I have assumed that the evolution of 
the fluctuations is described by the cold-dark-matter transfer function [32] where the 
present Hubble parameter is taken to be Hq = 50 km s _1 Mpc _1 . 

For power-law inflation, the best fit is given by n = 0.9 (bold curve in Fig. (lb)). 
However, n = 0.8 gives the best fit for natural inflation. From the galaxy data, 
there is not much difference between the best fits of these two models. One hopes 
to discriminate these models further when intermediate angle microwave background 
fluctuations [33] become more precise. Model 3, double inflation model is not partic- 
ularly attractive since it requires three parameters whereas the previous two models 
each required one less. For the choice of double inflation parameters advocated by 
Peter et al [34], there is not much advantage over the simpler models of power-law 
inflation and double inflation [2] (see Fig. (3b)). 



5.1 Disagreement with Grishchuk 

In a previous Erice meeting held in the fall of 1994, Grishchuk [35] stated that the 
tensor fluctuations arising from virtually all inflation models dominate the contribu- 
tion to the microwave background anisotropy observed by COBE [30]. I disagreed 
with this claim [36]. More recently, Deruelle and Mukhanov [37] have demonstrated 
in careful detail that Grishchuk was in error. Here I wish to clarify the issue further. 

The main point of disagreement arises during the computation of scalar pertur- 
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bations and not the tensor perturbations. Once again, the scalar perturbations are 
parameterized by the variable ((x) which is independent of time when the wavelength 
of the fluctuation exceeds the Hubble radius. Grishchuk's equation 

Cmfb = 7-7 (45) 
2n z 

on p. 225 of ref.[35] is incorrect. It should read 

Cmfb = (x - n 2 ^=) j (2n 2 ) (46) 

Setting X = 0, and then taking n — ► 0, one finds the standard result that 

Cmfb = -\ J ^ (47) 

approaches a constant which is not equal to zero as n — ► 0. (An exact solution of 
the perturbation equation demonstrates this point quite clearly; see, e.g., ref. [1].) 
Grishchuk erroneously concluded that (mfb = because his method of computation 
was too crude and he neglected the term described above. The scalar contribution 
to COBE's signal is indeed proportional to (mfb, an d for many inflation models, it 
dominates over the tensor component. 



6 Summary 



Hamilton- Jacobi theory for general relativity provides some deep insights into the 
structure of semiclassical superspace which now far exceeds investigations in homo- 
geneous models. Superspace describes an ensemble of evolving universes, and its 
complexity strains the imagination. However, the gradient expansion allows one to 
separate superspace into an infinite but discrete number of manageable pieces which 
are relatively easy to understand. 

The question of time choice in general relativity is a difficult one, particularly 
for the quantum theory [28]. For semiclassical problems of interest to observational 
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cosmology, one may construct a straightforward covariant formalism which treats all 
time choices on an equal footing. Different contours of integration in superspace 
correspond to different time foliations and they all yield the same answer for the 
generating functional provided that spatial gauge invariance is maintained. 

Although not quite perfect, reasonable fits of microwave background anisotropics 
and galaxy clustering may be obtained by power-law inflation with a spectral index 
of n = 0.9, or by natural inflation with a spectral index n = 0.8. If one wishes to pay 
the price for an additional parameter, double inflation is also adequate. 
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8 Figure Captions 



Figs, (la), (2a), (3a): Shown are the power spectra V^(k) for zeta, which describes 
the primordial scalar perturbations arising from inflation. Three plausible models of 
inflation are considered: (1) power-law inflation, (2) natural inflation and (3) double 
inflation. The first two models require two parameters: an arbitrary normalization 
factor and a spectral index n, where n — 1 is the slope of the spectrum in a log-log 
plot. Double inflation is a three parameter model. For each model, the normalization 
is fixed by large angle microwave background anisotropics. 

Figs, (lb), (2b), (3b): For the present epoch, the power spectra for the linear 
density perturbation Sp/p are shown for the same models of Fig. (a). The data 
points with error bars are the observed power spectrum derived from galaxy surveys. 
For power-law inflation, the best fit (bold curve) is obtained with a spectral index 
of n = 0.9, whereas n = 0.8 yields the best fit for natural inflation. (Gravitational 
waves are important for power-law inflation but not for natural inflation.) Using an 
additional parameter, double inflation also gives a reasonable fit. 
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Fig.(lb) 
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Fig.(2a) 
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Fig.(2b) 
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